arXiv:1509.07778v2 [math.AP] 14 Oct 2015 


Regularity of the velocity field for Euler vortex patch evolution 


Daniel Coutand 
Department of Mathematics, 
Heriot-Watt University 
Edinburgh, EH14 4AS UK 

email: d.coutand@ma.hw.ac.uk 


Steve Shkoller 

Department of Mathematics 
University of California 
Davis, CA 95616 USA 

email: shkoller@math.ucdavis.edu 


October 14, 2015 


Abstract. We consider the vortex patch problem for both the 2-D and 3-D incompressible Euler equations. 
In 2-D, we prove that for vortex patches with Sobolev-class contour regularity, fc ^ 4, the velocity 

field on both sides of the vortex patch boundary has regularity for all time. In 3-D, we establish existence 

of solutions to the vortex patch problem on a finite-time interval [0, T], and we simultaneously establish the 
j^k-0.5 j-ggyiarity of the two-dimensional vortex patch boundary, as well as the regularity of the velocity 
Helds on both sides of vortex patch boundary, for fc ^ 3. 
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1 Introduction 

1.1 The incompressible Euler equations 

Global existence for the Euler 2-D vortex patch problem was first established by Chemin 
Bertozzi & Constantin [2], and Serfati m-, see also the recent article by Bae & Kelliher [1]. Local 
existence for the 3-D vortex patch problem was proved by Gamblin & Saint Raymond [8]. 

We are interested in the regularity properties of the velocity field associated to the vortex patch 
evolution. In particular, we analyze the incompressible Euler equations on R", n = 2, 3, written as 

-I-Vp = 0 , (la) 

divu = 0 , (lb) 

where u(x, t) is the velocity vector field and p{x, t) is the pressure function, where the advection 
term denotes ■ 

1.2 The 2-D vortex patch problem 

Letting V-*- = (—we define the 2-D vorticity function uj{x,t) = V-^- ■ u{x,t) = m^,i —m^, 2 . 
The vorticity co is transported and satisfies 


iOt -f — 0 . 


( 2 ) 


Letting 'ij}{x,t) denote the stream function, given by m = we have that A'f/; = w, so that 

'tp{x, t) = ^ ^^2 log |a; — y\ijj{y)dy. Thanks to the Biot-Savart kernel K{x) = log |a;|, 

u{x,t) = r K{x - y)uj{y)dy . (3) 

Jr2 


Eor each time t e [0,oo), let D+(t) denote an open, simply-connected, and bounded subset of 
with boundary r(t) := dSl+(t) given by a closed curve which is diffeomorphic to the circle Let 
D“(t) denote D+(t) . The 2-D vortex patch problem consists of the following initial data for the 
Euler equations: 


UJo{x) 


1, a::GD+(0) 

0, xeD-(O). 


(4) 


The time-dependent open set fl'^it) is thus termed the vortex patch; the vortex patch boundary 
r{t) := dD+(t) moves with the velocity of the fluid, given by u{x,t) = K{x — y)dy. It follows 
that 

Vu{x,t) = I VK{x — y)dy. (5) 

Given an initial 2-D vortex patch boundary r(0) of Holder class it was established by 

Chemin [5] and Bertozzi & Constantin [5] that a unique solution exists for all time, that the 
contour regularity propagates, and that the gradient of the velocity remains bounded for all time. 
Their proof of contour regularity (in 2-D) can also be used to establish contour regularity 
(we provide a proof for the n-dimensional case, n = 2 or 3 in Section [5]), and we state one of their 
fundamental results as follows: Given an initial vortex patch boundary r(0) of class fc > 3, 

for all t e [0, oo), there exists a unique solution to the vortex patch problem, with non self-intersecting 
boundary r(<), and satisfying the following estimate: 


+ l|2(^ + !|V■^^(•J^)llL'»(R2) < F{t) , 


( 6 ) 
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where 

m, 


r(t) denotes an °-®-class parameterization of the vortex patch boundary 


\zU{t) 


inf 

Si¥=82 


\z{di,t) - z{e 2 ,t)\ 
\ 01 - 02 \ 


(7) 


and 0 < F{t) < co for any t < co. We see that (|6]) provides a strictly positive lower-bound on |z|*(t) 
which, in turn, provides a strictly positive lower bound for the metric \d0z{9)\ and ensures that T(t) 
does not self-intersect (see, for example, Majda & Bertozzi [in]). We identity with the interval 
[ 0 , 2 ^]. 


1.3 The 3-D vortex patch problem 

In three space dimensions, the 3-D vorticity w = curlrt is a vector field, and satisfies the vector 
equation 

ujt +^u} = , ( 8 ) 

where in components and for each i = 1,2,3, [V^w]* = X;j=i ■ 

Letting denote the vector stream function, given by u = — curl'0, we have that A'0 = w, 

and hence It follows that 

u{x,t) = IC{x - y)w{y)dy, (9) 

Jb2 

where IC{x) = is the Biot-Savart 3x3 matrix kernel. 

What type of vortex evolution in three space dimension is analogous to the 2-D vortex patch 
problem? The answer is as follows: we suppose that at time t = 0, D+(0) denotes an open bounded 
subset of which is diffeomorphic to the open unit ball B = {x e : |a:| < 1}. We then 

let r(0) = (5D+(0), and define D“(0) = D+(0) . We choose an initial divergence-free velocity field 
uo(x) = Uq Uq (a:)lf 2 -(o) such that the initial vorticity vector wq = curlitg G and 

satisfies 


, , ( curlu;i'(x), xgD+(0) 

curlew, ■ 

(10a) 

[wo • n(.,0)] = 0, 

(10b) 

[wo X n(.,0)] 7^ 0, 

(10c) 


where n(.,0) denotes the outward unit normal to dD+(0). The velocity uq is continuous across 
r(0) while the tangential components of ujq are discontinuous. The 3-D analogue of a 2-D vortex 
patch amounts to choosing uq in such a way that curlrtg = 0 on D“(0) and hence, necessarily, 
curlwQ • n(0) = 0 so that ojq is tangent to r(0). 

To explain this analogy, we first state the following existence theorem for the Euler equations 
m with initial data u{x,0) = uq{x). Gamblin & Saint Raymond [5] proved that whenever r(0) 
is a e (0,1), uq e 1 < p < go, and wq e 1 ^ q < 3 such that ujq has 

regularity in directions tangent to r(0), then there exists a unique solution u e L”(0, T ; IT^’“(IR.^)) n 
IT^’°°(0, T; LP(K^)) to ([1]). Eurthermore, letting r]{x,t) denote the Lagrangian flow of u, so that 

dtr]{x,t) = u{r]{x,t),t) for t>0, (Ha) 

r]{x, 0) = X , (11b) 

and for each t e (0,T], setting r(t) = r]{r{0),t), then r(t) is a closed surface of class and 
uj{t) e L'?(R^) such that oj{t) has regularity in directions tangent to T{t). 
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For each t G [0,T'], the Lagrangian flow is a diffeomorphism with Jacobian determinant 

det'Vr]{x,t) = 1. We set r2+(t) = ? 7 (n+( 0 ),<) and n~{t) = r]{Q~{0),t). Integrating the vorticity 
equation ([ 8 ]), we see that 

= Vr]{x,t) ■ ojoix) , (12) 

where in components, [Vr? • Wq]* = 

We will set the 3-D vortex patch problem inside of a periodic box. We let D denote a periodic box 
[—£, £\^ in with opposite sides of the box identified with one another, and with I taken sufficiently 
large so that D+(0) c D. Functions defined on D are 2£-periodic in each of the three coordinate 
directions, i.e., 

u{x + 2£ei) = u{x) 'ix e i = 1,2,3, 

were e\ = (1, 0 , 0 ), 62 = ( 0 , 1 , 0 ) and 63 = ( 0 , 0 , 1 ). 

The 3-D vortex patch problem has the following initial data: 

r( 0 ) is a closed surface diffeomorphic to , 

is an open set diffeomorphic to the unit ball in , 

D"( 0 ) = D-D+( 0 ), 

Uo{x) = + Uo_(x)la-(o) , 

div uq = 0 , 

ojQ = curl uq , 

/ , [ curlMn(a^): xgD+( 0 ) 

“"‘*>-(0, xen-(O) ' 

curlwQ • n(-, 0 ) = 0 on r( 0 ), 
curlwQ X n(-,0) 7 ^ 0 on r(0), 

uo{x)dx = 0 . 

Jn 

We then call r2’''(0) the initial vortex patch and r(0) the initial vortex patch boundary. The identity 
(fT^ shows that for each t e [0,T], uj(-,t) = 0 in D“(t) and that uj{-,t) ■ n{-,t) = 0 on r(<). In 
particular, if the initial vorticity is supported in a set which is diffeomorphic to a ball, then the 
vorticity stays supported in a set diffeomorphic to a ball for all time t G [0, T] for which the solution 
exists. In (fl^ . we could instead set f2“(0) = R^ — D+(0). 

1.4 Statement of the main result 

Because of the singular nature of VitT, it is difficult to establish regularity for higher-order derivatives 
of u with the formula ([5|) . By taking a different approach, however, we shall prove that the velocity 
field indeed enjoys higher-order Sobolev regularity on both sides of the vortex patch boundary. In 
particular, for the 2-D vortex patch problem defined in Section lOl we have the following 

Theorem 1 (Regularity of velocity field in 2-D). Given initial data and a global-in-time solution 
to the 2-D vortex patch problem satisfying 

■j—j— 7 — -I- ||z(-, t)||//fc-0.5(Sl) -I- II Vm(-, t)|| 2 ,<»(][j 2 ) ^ F{t) 

for t G [0, go) and fc > 4, the velocity field satisfies G iJ^(D+(t)) and u~{-,t) G Hlf^(Jl~{t)), 

and 

11'*^ (•j^)llff*^(n+(t)) + 11^ {n-{t))uB{o,R(t))) < G{t ), 


(13a) 

(13b) 

(13c) 

(13d) 

(13e) 

(13f) 

(13g) 

(13h) 

(13i) 

(13j) 
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where B{0,R{t)) is a ball centered at 0 with radius R{t) > 0 such that r(<) c B{0, R{t)), and 
G(t) > 0 is a function of F{t), defined in (0), with G(t) < oo for any t < co. 

Remark 1. Notice that both velocity vector fields and u~ gain a half-derivative of regularity with 
respect to the regularity of the vortex patch boundary T(t). This is very natural in Sobolev spaces 
, but requires us to locally extend our 1-D parameterization z{-,t) to a 2-D local diffeomorphism 
0+ (•, t) and 9~ (•, t) which also gains a half-derivative of regularity. This is accomplished by a specially 
chosen elliptic extension which we describe in Section [31 On the other hand, if we had assumed 
instead that the parameterization z{-,t) G then a standard local “graph” extension would 

have sufficed. More specifically, if z{-,t) is given locally by the graph {xi,h{xi)), then {xi,X 2 + h{xi)) 
provides a local extension to a diffeomorphism, but does not gain a half-derivative of regularity. 


Remark 2. Without any change to our proof, the initial data (jH) can be replaced by the more 
general initial data 


for any functions ujq 


UJo{x) 


Wq (x) , X G r2+(o) 
W 1 q{x) , X G ri“(0) 


G and ujq g ), A. 


Remark 3. In fact, Theorem |T] is true for fc > 3, but the proof requires one less regularization step 
for k ^ 4. 


Whereas Chemin |4] and Bertozzi & Constantin [2] have established regularity of the contour r(t) 
for the 2-D vortex patch problem, the regularity of the 3-D vortex patch boundary T(t) is limited 
to in the analysis of Gamblin & Saint Raymond [8]. As our final result, we simultaneously 
establish an existence theory in Sobolev spaces for the 3-D vortex patch problem, as well as the 
Sobolev-class regularity of the 2-D closed surface T{t) and the velocity fields and u~. 

Theorem 2 (Existence and regularity for the 3-D vortex patch boundary and velocity fields). For 
fc > 3, */r(0) is a closed surface of Sobolev class , and uq g F[^(il) with Uq g iJ*(D+(0)), 

Uq g and satisfying USl) . then there is a time T > 0 such that there exists a unique 

solution to the 3-D vortex patch problem, and for each t G [0,T], the vortex patch boundary T(t) is 
in \ u+{-,t) G H'^{n+{t)), and u-{-,t) G H^{n-{t)). 

Remark 4. The more general initial data (flUl) can replace (1131) in Theorem |3J 

Notation. We will denote the partial derivative by f,j for j = 1, 2, or 3. We will use the 
Einstein summation convention, wherein repeated indices are summed from 1 to n, with n equaling 
either 2 or 3. 


1.5 Outline of the paper 

In Section 121 we define the strong form of the two-phase elliptic problem that the two-dimensional 
stream function must satisfy, and we also define the associated variational formulation. In Section [31 
we define the local diffeormorphisms that we use to locally flatten the vortex patch boundary; these 
diffeomorphisms gain one-half derivative of interior regularity in F[^ spaces relative to the regularity 
of the vortex patch boundary. Section [4] is devoted to the Sobolev regularity theory of the fluid 
velocities u^{-,t) and u~{-,t) in the 2-D vortex patch problem. 

The 3-D vortex patch problem is studied in Section |S] After defining the two-phase-elliptic 
problem for the fluid velocity, we simultaneously prove existence of solutions and establish the 
regularity theory for both the vortex patch boundary r(t) and the velocity fields u~^{-,t) and u~{-,t); 
this is done in the Lagrangian framework. Einally, in Section |6l we establish the fundamental 
regularity estimates for the two-phase elliptic problem with Sobolev-class coefficients in n-dimensions 
(which arises in many applications, including the vortex patch problem). For completeness, we 
include a short appendix with some basic inequalities that are used in Section [31 
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2 A two-phase elliptic problem for the 2-D stream function 


The 2-D vortex patch problem has been previously studied using the evolution equation for the 
parameterization of the contour z{-,t) ([21 [5]); see also [5] for perturbations of circular patches and 
[B] for elliptical patches). We will take a different approach. 

While not necessary, it is convenient to introduce the stream function formulation of the problem. 
Let ^ — F~ on r(t), and let denote 

the outward unit normal to r(t), and T(-,t) denote the unit tangent vector to r(t). 

For each time t e [0, oo), the bounds ([1]) show that Vm(-, t) e L® (K^); thus, u{-, t) e and 

so the stream function G is a solution to the following two-phase elliptic problem for 

each fixed t e [0, oo): 


= in n(t)"'", (14a) 

Atp~{-,t) = 0 in n(t)~, (14b) 

[V'(-,1)1 = 0 on r(t), (14c) 

=0 on T{t). (14d) 


The fact that ^/’(■,t) G means that the interface jump condition ([MU') holds in 7J°-®(r(t)). 

For each time t G [0, oo), (imi has the following weak formulation: 

I ■V<j)dx+ \ -Vcjidx = - \ <j)dx y(l)e . (15) 

Ja+(t) Jn-(t) Jn+(t) 

From the bounds the stream-function satisfies 

\\'^{-,t)\\H^{BiO,R{t)) < F{t) , (16) 

where 5(0, R{t)) is a ball centered at 0 with radius R{t) > 0 such that r(t) cz 5(0, R{t)). 


3 Locally flattening the boundary r{t) 

We construct local diffeomorphisms in small neighborhoods of r(t) which locally “flatten” the vortex 
patch boundary, and which gain one-half derivative of regularity in the interior with respect to 
the regularity of the parameterization There are other methods to construct regularizing 

diffeomorphisms (see, for example, [21171 UDi but the method we present appears quite natural for 
arbitrary geometries. 

Let 5+ = {x G : |x| < 1} denote the open unit ball in with boundary = {x G R^ : 
|x| = 1}, the unit circle. For each t G [0, oo), we solve the following elliptic equation for Z(r, 9, t): 


A^Z+ = 0 in 5+ , 

(17a) 

Z+ = z on 

(17b) 

dZ+ dz-L , 

dr ee ^ • 

(17c) 

The unique solution Z'^{r,9,t) to (1171) satisfies the estimate 

■,t)\\Hk(D + ) < C' ^(’jA) j 

(18) 


and we are considering integers A: > 4. The boundary conditions (Hlb ,c) show that 

detVZ+(l,6i,t) = \dez{e,t)\^ . 


6 






D. Coutand and S. Shkoller 


Regularity of the velocity for Euler vortex patch evolution 


From the definition 0 of \z\^{t) and its lower-bound given by dH), it is proven in [10] that there 
exists a function a{t) > 0 such that a{t) sg minggsi \dgz{9,t)\‘^. Hence, det VZ+(1,0,t) > a{t) > 0. 
This shows that is locally injective around each point on S^. 

Next, we define D~ = {x e : 1 < |x| < R{t)}, where R{t) > 0 is chosen sufficiently large so 
that the ball B{0,R{t)) contains r(t). We let Z~{r,9,t) solve 


A^Z- = 0 

Z- = z 

z- = Id 
dZ- _ dz^ 
dr do 
dZ- 


in D~ , 
on , 

on [r = R{t )}, 
on , 

on [r = R{t )}, 


(19a) 

(19b) 

(19c) 

(19d) 

(19e) 


where denotes the unit basis vector (cos 0, sin 9). Again, we see that the unique solution Z {r,9, t) 
to (fTOjl satisfies the estimate 


(•ri^)llH'=(D-) < C'l|2(’i^)llff'“-o-5(si) ■ ( 20 ) 

We define the map Z = Z'^l-^+ Z~lu-. Due to the boundary conditions (fTTb.cl and (IT^.d') and 
the Sobolev embedding theorem, the map (r, d) Z{r,9,i) is and for any point 9 e there 
exists a ball B{9,e{t)) a R^, centered at 9 with radius e{t) > 0 taken sufficiently small, such that 
Z{-,-,t) is injective on B{9,e{t)). 

Next, we show that for e > 0 sufficiently small, the image Z+(l — e,9,t) is contained in H+(t), 
and similarly, that the image Z~{r,9,t) is contained in f2“(t). To that end, let 0o(^) denote the 
point in [0, 27r] at which the maximum value of z(9,t) ■ ei occurs. We assume that the tangent 
vector dez{9o{t),t) = f3{t) ei for some /3(t) > 0 (for, otherwise, we can reverse the orientation of the 

I dZ'^ 

parameterization). Hence, dgz-^{9Q{t),t) = j3{t)e2- This shows that -^^{l,9Q{t),t) > 0, which in 
turn implies that Z^(l — e,9o(t),t) < Z^(l,9o(t),t) which proves that, for e > 0 sufficiently small, 
for all r G [1 - e, 1) and 9 e [Ooit) - e, 9o{t) + e], 

Z"^{r,9,t) ■ 62 < z{9o{t),t) ■ 62 . 


Therefore, maps a local neighborhood of 9o{t) (in D'^) into H+(t). Since Z^ is locally injective 
around §i, this means that the image of any Z’+(l — e, ■,t) (for e > 0 small enough) stays in H+(t), 
otherwise it would intersect r(t), which we shall next prove that this cannot occur. Similarly, the 
image of any Z~(l -I- e, ■,t) stays in 

We next prove that for e > 0 sufficiently small, 

Z+{1 - 6, 9, t) n r{t) = 0 V6» G 

Since det VZ+(1,0,t) > a{t) > 0 for all 0 G by the inverse function theorem, there exists a 
small ball B{9,TZ{9)) c centered at 0 G with radius TZiO) > 0, such that Z+(-,-,t) is a C^- 
diffeomorphism between D’*' n B(9,71(9)) and Z~^(D''' n B(9,7Z(9)),t), as well as a homeomorphism 
between D+ n B(9,7Z(9)) and Z^(D+ n B(9,7Z(9)),t). Since the compact set is covered by 
IJgggi i3(0,7^(0)), we can extract a finite subcover \J^^^B(9i,7Z(9i)), where 9i, i = are 

points in . 

Let = {a; G R^ : 1 — e ^ |a;| < 1} denote an annulus. We choose e > 0 small enough 
so that 21*^ c (_jf^^B(9i,7i(9i)). With (r,9) e A'^ fixed, we choose i e such that 

(r,9) G B(9i,7Z(9i)) n Z1+. Since Z+(-,-,t) is diffeomorphism between I?+ n B(9i,7Z(9i)) and 
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n B{9i,TZ{9i)),t), then Z'^{r,9,t) e Z+(Z?+ n B{9i,TZ{9i)),t). Furthermore, as is 
an homeomorphism between D+ n B{9i,7i{9i)) and Z+(Z3+ n B{9i,TZ{9i)),t), then Z+(r,0,t) ^ 
Z+(d[-D+ n B{9i,TZ{9i))],t), which implies that 9, t) ^ z{\9i — TZ{9i),9i + TZ{9i)],t) cz r(t). 

In summary, we have shown that for {r,9) e B{9i,TZ{9i)) n _D+, 


Z^(r,9,t) is in the interior of Z^{D+ n B{9i,TZ{9i)),t) 


with 

diameter sS 2||VZ+||ioc(I5+)7^(6»i). 

From the positive lower bound ([5]) on the function |z|*(i) in (O, there exists eo > 0 such that 
for any x e r(t), B{x, cq) does not contain any point of r(t); therefore, choosing the radius 

'R-{9) such that 

2||VZ+||j;^Qo(£)+)72.(0i) < eo , 

(and increasing N if necessary) we have that Z^{D^ n B{9i,TZ{9i)),t) does not contain any point 
of r(t), which shows that Z'^{r,9,t) ^ r(t) as desired. A similar argument shows that for {r,9) e 
B{9i,TZ{9i)) n D~, Z~{r,9,t) is contained in n~(t). 

Thus, for each 9i e i e {1,...A^}, let Ui{t) = B{9i,TZ{9i)) cz and let Vi{t) = Z{Ui{t),t). 
The map Z is then a diffeomorphism oiUi{t) onto Vi{t), and due to the estimates (ITSl) and (l20l) . 

Z'-(-, •,t) : D- n Ui{t) n Vi{t) is an diffeomorphism. 


Next, we flatten the boundary of h{i{t) n For each i e {1,..., A^}, h(i{t) n is a graph given 
by {xi,hi{xi,t)) where each hi{-,t) is We define the C® local diffeomorphisms 'd^{t){xi,X 2 ) = 
{xi,X 2 ± hi{xi,t)) with det = 1, and we set 

BX = n D^) and RJ = [^?+(<)]-i(W,(t) n §') . 


The set B^ cz {x 2 =0} is a flat boundary. 

Finally, we define 9^{t) = Z-{t) o Then 

9f{t) : B- fl- n Vi{t) is an diffeomorphism, 

and thanks to (HU), dsni), and ®, for each i e {1,..., N}, 

1 


det V6i±(t) 


L®(B+) 


+ ll^i (^)IIb'=(b+) < 'PiF(t )), 


( 21 ) 


( 22 ) 


where 'P{F{t)) denotes a generic polynomial function of F{t). Furthermore, if we set 9i{t) = 
9^{t)l-^ + 9~{t)lgi , then each 9i{t) e where B = u i?_ u Bq. 


4 Regularity of the velocity field for 2-D vortex patches: 
Proof of Theorem [1] 

We first use the weak formulation CSl) to build regularity of the stream function '0-. Interior 
regularity of 'k- on sets away from the patch boundary r(t) is classical, so we focus our attention 
on regularity of near r(t). We will use the change-of-variables 9i{t) given in (I^Tl) . 

Step 1. The elliptic problem for ^k- set on B+. The weak formulation (IT5|) can be written as 
[ V0;"'"(-,t) • V0da; + [ ■ Vcjidx = — [ (j)dx 
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for all test functions </> e i?Q(V(t)) and each 

With the collection of diffeomorphisms {Oi}fLi given in (|2T]) for each t e [0, oo), we define 


Af = [Vd±(t)]~i and J±(t) = det Vd±(t) 


and set 


At = ■ 

It follows from (l22l) . and ([61) that for all t e [0,oo), there exists a function 0 < Xi{t) such that 

w'^ Afix) w ^ Xi{t)\w\'^ yw e M? , x e B^.. 


(23) 


To establish (1^ . we drop the i subscript (and superscript), and let w+ = J+'^A+ w. The left-hand 
side of (1^ is simply |w+P, and w = therefore, 




|| J +'/" V 0±||2 




L>»(B + ) 


SO that A(t) = which has a strictly positive lower bound since X(t) ^ = 


J+^^^Vd-||^co(B+) < 'PiFit)) by (l2^ . Additionally, from 


l|Al+||_H-fc-i(_B+) < CV{F{t)) . 


(24) 


We set 


^f± = i/,±o 6 », $ = (/)o 6 ». 


Since ij) e HQ{V{t)) and each 9i{t) e '^^(R), it follows that $ e Hq{B), and can thus be used as a 
test function. By another application of the change-of-variables formula, we then have that 






+ k (•,!) dx + 


( A*'A'^~,k (•, t) dx = -\ $ J+ dx V$ e Hl{B). 

JB- JB+ 


(25) 


Step 2. regularity for -0+ and -0 . We set fc = 4 so that 0- e H*{B+) and first establish 
that each ijj- is H^. We let (Cilfli denote a smooth partition-of-unity, subordinate to the open 
cover in particular, 0 ^ Ci ^ 1 in denote a smooth cut-off function, C* = 

and let = Q o0^{t). We define the horizontal convolution operator as follows: for e > 0 sufficiently 
small, 

AeE = Pe{xi - yi)F{yi,X2)dyi , 

Jr"-i 

where Pe{xi) = e“^p(xi/e), and p is the standard mollifier on K. We again drop the i subscript, and 
substitute 

$ = e H^{B ), 4- = 

into (HH). Since differentiation commutes with convolution, we have that 


The variational formulation (1^51) then takes the following form: 


j± 4-12* = - f d? j+ e^idfie 'I') dx , 
Jb. 


(26) 
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where 

I± = f ),11 A2(e2^±),,n dx, 

JB± 

= -2 f mj dx. 

JB+ 

Next, we see that 

J± = r A,(e2^±)„n dx+ f {{A,,A’^^}ie^^),kii ) A,(e2^±)„n dx 

JB+ JB+ 


+ J A, [24^1 (e^vl;),,, +A’^^ai A,(C2vI/±),^.ii dx , 

^-V-" 

where 

{A„ 4 ^}(e^M/±),fen = A,iA'^^{e^±),kii ) - A,(e'vi,±),,ii (27) 

denotes the commutator of the horizontal convolution operator and multiplication by A^. Using 
the lower-bound (l23t . we see that 

A(t)||d?A,V(fM/±)||i.(5^) (28) 

We let 0 < (5 « 1; we will make use of the Cauchy-Young inequality ab ^ SX{t)a‘^ + 

a,b^ ^■ 

Using Holder’s inequality together with the Sobolev inequality ||/||lp(b+) < C\\f\\H^(B+) for all 
/ e H^{B+) and all p e [1, 00 ), we have that 

< ^^11-4+ ,11 IIl2(B+) ■ 

Thanks to (l24l) and (fT6l) . we then infer that 

|Xi±|^7^(U(t))||A,(f^k±)„n|U=(B,) 

6XmUe^^),jn Wh^BA + 

^ dA(t)||A,(e^vi;±)„n \\h(BA + (1 + {SX{t))-^)V{F{t)), 

where we continue to use V{F{t)) to denote a generic polynomial function of F(t). A similar estimate 
can be established for the integral , which provides us with 

|Ii±| -b |I±I [1 + {5X{t))-^\P{F{t)) + 5X{t)\\A,{iH^)Bi Vh.^ba • (29) 

Also, the integral on the right-hand side of (1^51) has the same upper bound. 

It remains to establish such an upper bound for |Xib|. We set then, 

[A,:(AI+-^ 5^) - A,p^](a;i,a;2) = [ pe(a;i - yi)[^+^(yi,X 2 ) - (xi,a;2)]5fe (j/i, X2)dyi. 

J xi —e 

(30) 
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From Morrey’s inequality, for all yi e B{xi, e), 

\A'^{xi,X 2 ) - A'^{yi,X 2 )\ ^ Ce sup \A'^,i{yi,X 2 )\ ^ Ce\\A\\H^(^B+) ■ (31) 

yi€(xi-e,xi+e) 

Substituting m into (1501) and using Young’s inequality for convolution, together with (|001) . we see 
that 

||{A„4^Kf ^ ^ Cer{F{t))\\V{e'f^),ii ||l=(b±) , 

so that 


\Iu\^CeV{F{t))Me'f^)ai\\UB,)- 


We choose e sufficiently small so that CeV{F{t)) < X{t)/2. By choosing <5 > 0 sufficiently small, we 
obtain from dUl), (EHl) and (1001) the estimate 



\^dx^ [1 + X{t)-^]V{F{t)). 


Passing to the limit as e ^ 0, we find that 



fdx^ [1 + X{t)-^]V{F{t)). 


(32) 


From (ITOh .bl. we have the following identity holding at any point of the interior of B+: 

,221 = 2^+^4'-,2ii +2A ^,i 'I '*,21 +A }^,1 'k*,ii 

+ Al^i ^±,22 +A^,ji +A!^,j ^±,fei . 


The lower-bound 


shows that A^ > X{t); hence from (1501) . (1051) . and (ITOl) . 



P dx ^ [A(t)-i + X{t)-^]r{F{t)). 


(33) 


Then, since 


- Af^-,222 = 2 Af ' i /-,221 +^+ «'*,211 + 2 A ‘+,2 ^'*,21 + A +,2 ^'*,11 
+ Af ,2 ^'-,22 +^±",22 '^-,k + A^±,J «'-, fc2 , 

we use (|33l) to conclude that 

f e '|^±,222 I" dx ^ [X{t)-^ + X{t)-^]r{F{t)). 

JB+ 


(34) 


Given the interior estimates, we sum ([32|) . (l33l) . and (l34|) over our finite cover index i = 1, 
and find that 

l|■*/’^(^^)ll^^-3(Q+(t) + \\ip (•>^)llff3(n-(t)nB(o,fl(t))) ^ ^(^(^))) (35) 

where we have used the fact that A(t)“^ sS V{F{t)). Then since u- = V-*-'0-, (ISSl) shows that 


ll^''"(^^)llff 2 (n+(t) + 11 ^ (’:^)llff 2 (Q-p)nB(o,fl(t))) ^ ■ 


(36) 
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Note that the estimate (IMI) has been obtained for the case that r{t) is of Sobolev class H^-^) so 
that we can indeed build further regularity for u-. 

Step 3. regularity for and u~. We will now use estimate (IMl) to build the regularity 
for u'^ and u~. On T(t), we let VrW denote the directional derivative of u in the direction t and 
similarly, we let denote the directional derivative of u in the direction n; for example, in 

components = u^,j ■ We make use of the following identities on r(t): 

divM = VrU • T + VnU- ' u , (37a) 

curln = VrU • n — 'S/nU ■ t . (37b) 

Since is continuous across T{t), it follows that 

[VnM • r] = -[curlul = - 1 . 

Then, using (I57h'). and the identity 

|[V„n] = |[V„M • T|r + |[V„M • n|n, 


we see that 

|[V„u| = -T. 

From (15^ , the velocity field u is a solution to the following two-phase elliptic problem: 

Au- =0 in , 

[m] = 0 on r(t), 

[V„m] = -r on r{t) 

with variational form given by 

r : Vwdx + r -.Vwdx = f T{-,t) ■ w dS{t) Vw e R^), (38) 

Jn+(t) Ja-(t) Jr(t) 

where A : B = for any 2x2 matrices A and B. 

Again dropping the subscript i, we write (1551) locally as 


I Vu''"(-, t) : Vw dx + I Vu (•, f) : Vw da; = I T{-,t) ■ w dS{t). (39) 

Jv(t)nn+(t) Jvft)nD-ft) Jv(t)nr(t) 


■'V(t)nQ (t) 


V(t)nr(t) 


for all w e idp(V(t);R^). We set U = uo9 and W = w o 0. By the change-of-variables formula, (155)) 
becomes 


f A+m+i-A) ■W,jdx+ ^ At^U-{-,t) -Wjdx = - \ 9,i-WdS yW e H^{B-R'^). (40) 

J JB — '-'^0 

We then substitute 

W = e^ldtieu) e H^,{B), 

into OMl). By repeating the identical argument of Step 2 above, we find that 

^)llfi-3(n+(t) + 11'“ (■>^)llfa3(n-(t)n_B(o,iJ(t))) ^ ■ (41) 

Step 4. H'^ regularity for and u~ . We continue to assume that fc = 4 so that the boundary 
r(t) is of Sobolev class and our change-of-variables 9^{t) e H^{B+). We will now show that 
u'^ and u~ have H'^ regularity. 
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To do so, we let the test function W = in (l40l) . By a slight modification of Step 

3, we find that 

+ ||m (■,i')HH4(n-(t)nB(o,Ji(t))) ^ • ( 42 ) 

There are new types of integrals that arise in establishing the regularity; namely, integrals that 
have highest-order derivatives on both C7- and 0-. 

One of these integrals is analogous to one of the integrals in Ii)r defined in Step 1 and is written 
as 

J± = r A, ieu^).k] dx. 

^ B+ 

We estimate the integral \J-\ using an L'^-L'^-L? Holder’s inequality: 

\J-\ SS ||^+ ,111 ||L2(B+)||(C^C4)-,fe ||l<»(B+) ||Ae(C^C/)-,jiii > 

which, with the Sobolev embedding of H'^{B±) into L'^{B±), shows that 

\J-\ < C'||^+ ,111 ||L2(B+)l|t^*,fc \\h‘^{B+) ||Ae(5^17)-,jiii \\l‘^(B+) ■ 

Using the estimate (l 22 l) with fc = 4 together with the previous lower-order estimate (l4T]) of it- in 
we obtain that 

\j±\^r{Fm\ueu^),mi iil^(b±), 

which is just a linear term in ||Ae(^^[/-),jiii \\l^(b+)j easily controlled by the energy integral 

Xl^ = X{t)\ |A,V(e"C7)±,iiipdx, 

JB+ 

analogous to the term Iiy in Step 2 above. 

Other integral terms set on B+ of this type arise and can be treated similarly. There is one slight 
variation: the boundary integral term 

Ie=\ dfA,{eO,i) ■ dlAeieU) dS, 

JBo 


for which we simply notice that 

\Xg \ ^ \\dlA^{^'^6,i)\\H-o.^(Bo)\\dlA^{^'^U)\\H0.5(Bo) 

< C\9b 11/1-3.5(Bp)||diAe(^^C/)||/i-i(B^) , 

where we have used the properties of the convolution operator for the first norm on the right-hand 
side, and the trace theorem for the second norm. This then provides us with 

|I,| ^iP(F(t))||A,(e2^±),iii Whhb^), 

which is a linear term controlled in a similar manner as 4 . 

Step 5. regularity for it"*" and u~. Letting W = in (HOI) and 

repeating Step 3, concludes the proof. □ 
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5 Existence and regularity of the 3-D vortex patch boundary 
r{t) and u±: Proof of Theorem [2] 

5.1 The two-phase elliptic problem for velocity 

As defined in Section ITT^ the vortex patch boundary r(t) is a closed 2-D surface which is diffeomor- 
phic to the unit sphere and that is an open subset of such that = T{t), and 

n~(t) = — n+(0). We denote by 12 in what follows, and we set 12- = 12-(0). 

We let and T 2 {-,t) denote an orthonormal basis of the tangent plane to each point of T{t), 

so that (ti, T 2 , n) is a direct orthonormal frame of We let Vt„u {a = 1,2) denote the directional 
derivative of u in the direction Tq, and similarly, we let denote the directional derivative of u 
in the direction n; for example, in components = u'‘,j (ray. We make use of the following 

identities on r(t): 


divit = • Ta + 'S/nU ' u , (43a) 

curlu = (V1-2W • n — VnU ■ T 2 )ti — • n — VnU • ti)t 2 , (43b) 


where we have used the fact that curlM+ • n = 0 n r{t) by (ICTrl. Since u{-,t) is continuous across 
r(t), it follows that 


|[V„M • Ti] = [curlw • T 2 \ = curlM+ • T 2 , (44) 

— [V„M • T 2 ] = [curlw • Ti] = curlM+ • ti . (45) 

Then, using (|43h l. and the identity 

[V„m] = |[V„M-Ta]T„+|[V„w-n]n, (46) 

we see that 

|[V„u] = curlw^ ■ T 2 Ti — curlrt^ • ti T 2 . (47) 

The velocity field u = rt-lQ-p) is a weak solution to the following two-phase elliptic 

problem: 


-Am+ 

= curl curin'*' 

in 

n{t )+, 

(48a) 

Au~ 

= 0 

in 


(48b) 


= 0 

on 


(48c) 

[V„n] 

= curln"*" • T 2 Ti — curin'*" • ti T 2 

on 


(48d) 


with variational (or weak) form given as follows: For all vector test-functions w e i7^(f2) given by 


■.'SI'w dx + Vu {■,t):Vwdx= curl(•, t) • curl wda: 

Jn+(t) Jn+(t) 

-I- [n X curln+] • wdS{t) 

Jr(t) 

-f I [curl(•, t) • T2 Ti — curlM^(-,t) • Ti T2] • i(;d 5 '(t), ( 49 ) 

Jr(t) 

where A: B = A'^Bj for any 3 x 3 matrices A and B. Next we notice that 

n X curlrt^ = n x [curlw^ • ti ti -f curin'*" • T2 r2] = curln^ ■ ti T2 — curln^ • 12 ti , (50) 
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so that the boundary integral terms of ( 1 ^ ) cancel each other, and we are left with 

: Vw dx + Vu~ {-jt) : Vw dx = curl u"*"(•, t) • curl re dx Vre e . 

Jn+(t) Jn-(t) Jn+(i) 

(51) 

We let r]{x,t) denote the Lagrangian flow of u, as defined in (fTTl) . We set v- = u- o rj and we 
define A{x, t) = [\/r]{x, Then, letting (j) = w orj, (ISTI) can be written as 

r dx + f dx = r [curl u"*"] 077 -Icurhf/) 077 “^)] o ?7 detVrydx (52) 

Jq+ dxj dxk Jq- oxj dxk Jn+ 

for all 4> e id^(f 2 ), where 

For solutions to the Euler equations (CD, divM = 0 so that detVry = 1 , but the general form (15^ 
will be necessary for our fixed-point scheme. 

5.2 The fixed-point procedure for existence of solutions to the vortex 
patch problem 

In Section [BJ we will establish the fundamental elliptic regularity results for a Lagrangian variational 
formulation as in (ED- Using that regularity theory, we now prove the existence and regularity of 
solutions to the 3-D vortex patch problem; our solutions have smooth Sobolev regularity on both 
sides of the vortex patch boundary r(t) and are globally in iJ^(D). 

5.2.1 The functional frame-work 

We remind the reader that we use D to denote a periodic box [—£, in with opposite sides of 
the box identified with one another, and with £ taken sufficiently large so that D+(0) c D. Functions 
defined on fi are 2 Uperiodic in each of the three coordinate directions, i.e., 

u(x + 2£ei) = u{x) Mx e i = 1,2,3, 

were ei = (1,0,0), 62 = (0,1,0) and 63 = (0,0,1). Functions in iL^(D) satisfy periodic boundary 
conditions, and iL^(f 2 ) can be identified with iJ^(T^). 

Given T > 0 and M > 0 assumed fixed, we work in the Lagrangian framework and define the 
bounded closed convex and nonempy set 

= {u e L"(0, T; H\n)) n T; iL'=(D±)); \\vh.^o,T-mm + ll^^±llL=(o.T;i7'=(o±)) < M} (53) 
for integers A: > 3. For any 7 ; e V^, we define the Lagrangian flow 

ri{x,t) = X + v{x,s)ds, (54) 

Jo 

which therefore, from (IBBll . satisfies rj e ‘^°(0, T; iL^(D)) n “^“(O, T; iL^(D-)). Note, also, that since 
the vortex patch boundary is transported by the fluid velocity, we have that 

r(t) = 7?(r,t). 

Hence, the regularity of the velocity field in £1'^ provides us with the regularity of rj in D"*"; the trace 
theorem then provides the regularity of 77 on L, and this in turn provides the regularity of the vortex 
patch boundary T{t). 
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Since is a periodic box, and hence convex, any two distinct points x and y in can be connected 
by the straight-line segment {x,y); therefore, by splitting the segment {x,y) into a finite union of 
subsegments {xi,Xi+i), we can assume that each subsegment {xi,Xi+i) is contained in either or 
n~. It follows from that 

r]{x, t) - r]{y, t) =x - y + v{x, s) - v{y, s) ds 

Jo 

=x — y + v{xi,s) — v{xk,s) ds 

Jo 

K-l 

=a;-y+ V s) - s) ds , 

^=l Jo 

which therefore shows by the fundamental theorem of calculus, that since each (xi,Xi+i) is either 
contained in or 0“, that 


n—1 

\r]ix,t) -r]{y,t) - x - y \ ^ C Y \xi - Xi+i\ || Vn(-, s)||i»(o+) || Vn(-, s)||l<»(q-) ds , 

i=i JO 

and from the Sobolev embedding theorem, 

ra-l pt 

\r]{x,t) -r]{y,t) - a; - j/| ^ C ^ \xi - x^+i\ || Vn(-, s)||/f 2 (Q+) || Vn(-, s)||/f 2 (n-) ds . 

i=i Jo 

From the definitions (155)) and (1551) . it follows that 

n—1 

|? 7 (x, t) - T]{y, t) - X - y\ ^ ^ \xi - \2ViM ^ 2-\/TMC\x - y\ . (55) 

2=1 


We now choose T such that 


so that for any x and y in 17, 


0 < T ^ 


1 


167 ^ 2(72 ’ 


(56) 


\-q{x,t) - y{y,t)\ ^ -|a: - y \, 


which establishes the injectivity of ry in 17. Furthermore, since 


|V77(x, s) — Id| ^ 


Vu(x,s) ds ^ ||Vn+(-,s)||ia,(a+)-t I|Vn (•,s)||z,co(o-) 

Jo Jo 


ds 




2CVtM , 


(57) 

due to the continuity of the determinant at Id in R®, we can choose T > 0 small enough, so that for 
all X e 17 and 0 ^ ^ T, 

^ > detV77(x, 7) i , (58) 

which shows, with the previously established injectivity, that ?7(-,7) is an diffeomorphism from 
17- onto the image ?7(f7-,7), and a homeomorphism from 17 onto 77(17,7). Finally, by choosing T 
sufficiently small we can ensure the strict positivity of the coeffiicient matrix A', for all 7 G [0,r], 


vY Af{x,t) w ^ ^ ■ 


G 17 . 


(59) 
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5.2.2 The fixed-point procedure 

We define the Lagrangian curl operator curl^ as follows: if t) is an Eulerian vector, and v = uorj, 
then we curl^ v = [curlu] o r] where for any differential vector field F, and for i = 1, 2, 3, 

[curl^ F]i = Eijk ’ (60) 

where denotes the permutation symbol, so that £ijk = 1 for even permutations, Eiji^ = — 1 for 
odd permutations, and £ijk = 0 otherwise. We will employ a fixed-point procedure on the variational 
equation (15^ . which we write as 


f aFjI 


f 



dxk 

Ja- 

OXk 


[curlu"*"] 

Jn+ 


o T] ■ curl^ (f) det Vridx 


(61) 


for all (j) e H^{n). From (I12|) . 


cmlu o rj = Vrj ■ ojQ , wq = curlitJlQ+. (62) 

Since divwo = 0, using the formula (1621) . we see that 

I curlw"''o ?7 da; = ( 77 (curluj • n(-, 0)) dS'(O) = 0 , (63) 

Jn+ Jr 

where the last equality follows from ([TSh). 

Now, given v in our convex set and letting ry denote the homeomorphism defined in (I54|) . we 
define 

C{v){x,t) = Vri{x,t) ■ uiq{x) in fl. (64) 

Notice that for any a; e E, the trace on E of C(v)(x, t) ■n{ri{x, t),t) (the trace taken from from 11+) 
is zero, and is thus equal to the trace of of C{v){x,t) ■ n{ri{x,t),t) evaluated from 11“. To see this, 
we use an important geometric property of the inverse deformation matrix A(x,t) = [V 77 (a;, t)]“^; 
namely, if N{x) := n(0,a:) denotes the outward unit normal to dH+ and if n{r](x,t),t) denotes the 
outward unit normal to dH+(t), then 


ni{r]{x,t),t) 


A'tNk 

\A^N\ ■ 


Hence, it follows that 


C{v) ■ nor] =C{vy 




\A^N\ \ATN\ 




(65) 


where we have again used (|T^1 for the last equality. 

Furthermore, the same computation as in (1631) shows that 


r C(v) dx = ( r] (curluj • N) dS{0) = 0 . (66) 

Jn+ Jr 

Now, for each time t e [0, T], we construct a solution v(-, t) to the following variational problem; 

I C(n) • curl,, 0 det VTyda; Vt/) e iJ^(ll). (67) 

Jn+ 


1 


1 

r 


Jn+ ^^3 

dxk 

Jn- dXj 

dxk Jr 
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From (l5^ and the Lax-Milgram theorem, there exists a unique periodic solution v{-,t) e for 

each fixed te [0, T], satisfying 

I V dx = 0 , (68) 

Jn 

Furthermore, since C{v) e k ^ 2, we may integration-by-parts on the right-hand side of (I67|) . 

We use the fact that the cofactor matrix a{x, t), defined by a = det Vr] A, satisfies the Piola identity 
= 0 for j = 1, 2, 3. Thus, we see that (l67l) can be written as follows: for all ij) £ 


1 




Jn+ 

dxk 

Jn- 

dxk Jf 


n+ 


-b 


I 

Jdn+ 


C{v) X {jN)(l)dS{Q), 


which is the variational form of the general elliptic system m studies in Section with forcing 
functions 


f_ = 0 , f+ = curl,, C{v) det Vr ], and g = C{v) x (a^N ), 
for which our regularity result Theorem 0] applies. We therefore have that (for fc > 2) 


||'(;+||^fc+i(Q+)-b ||n ||_f/fc+i(o-) 

<C[||/+|| + Il9llff'=“0-5(r) + "P(M+llH>'(n±)) (^ll/+llL2(n±) -b ||g||/i--o.5(r)^ j , (69) 

where 7^ is a polynomial function and the constant C depends on 
From (l57l) and for 

Vt M ^eo, (70) 

with 0 < eo « 1 denoting a sufficiently small constant (which is independent of M), that for any 

ll^llff‘^+i(n+) ^ C'|n|. (71) 

Since from the definition (IMl) . 

l|C(n)||//fc(o+) ^ ||wo||ii-fc+i(n+)(l + CVtM) , (72) 

we then infer from (l72]l . (fTTll and (j69l) that 

||h'''||j^fc+i(Q+)-b llti ll/ffc+qn-) < C'|^C|n|||Mo||fffc+i(n+)(l + C'eo)(l + '^(l^l))j • 


Therefore, 

which thanks to dZOl) shows that 

ll^''"llL2(0,T;ff‘^ + i(a+)) + l|f' llL2(0,T;H'= + i(n-)) 
This inequality then proves that n e for 




2c[c|r!|||?ro||ff'^+i(ne)(l + Ceo)(l + iP(|r!|))] Vt, 


^ 2c[c|f2|||nollff'=+qn^)(l + Ceo)(l + iP(|f2|))] ^ . 


= 2c[c|f2|||no||ff^+i(n+)(l + C^eo)(l + T’d^^D) 


eo ■ 


Moreover, it is easy to check that the map 0 : u i—> i; is sequentially weakly lower semi- 
continuous; that is, if Vj ^ n in the weak topology of the norm defining the closed convex set V^, 
then Qvj Qv. Therefore, by Schauder’s second fixed-point theorem (see [13], page 452), which is 
itself a corollary of Tyhonov’s fixed-point theorem, we then have that 0 has a fixed point in V^. 
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5.2.3 The fixed point is a solution to the Euler equations 

We now explain why this fixed point, n = u, is indeed a solution of the Euler equations with initial 
data Uq, and hence a solution to the 3-D vortex patch boundary. At a fixed point v = v, dSZl) 
becomes the following variational problem: 


I 




dv'^ dfj) 


n+ dxj dxk 


dx+ { 

Jf 


dv d(j) 


o- dxj dxk 


dx = 


f C{v) 
Jo+ 


curl^ (/) det Vry dx (73) 


where the operator curl^ is defined (l60l) . We define the following Eulerian quantities associated to 
our Lagrangian velocity v and test function (j): 

u = V o ri~^ , = C{v) o , and w = <f>or]~^. 

The change-of-variables theorem shows that (1751) can be written a^ 


I 


?7(f2+ ,i) 


• Vw dy + 


I 


ri{Q ,t) 


Vm • Vw dy = 


I 


r)(n+.t) 


£ • curl w dy. 


(74) 


Our goal is to show that divit = 0 and that curlu = £. To do so, we use integration-by-parts on 
the left-hand side of (1751) : we see that 


Vu • Vw dy = — \ Au ■ w dy + • w dS{t) 

= curl curl u ■ w dy — \ V div u ■ w dy + • w dS{t) 

Jr}(0,.t) *Jr}(0,.t) JriiT.t) 


I 




curl u ■ curl w dy + 


I 




div u ■ div w dy 


+ (|[V„m| -I- |[n X curln] — |[divM|n-) • wdS{t). 

JvlT.t) 


h{r,t) 

The identities (l46l) and (l50l) show that for u e so that |[m] = 0 on we have that 

-I- |[n X curln] — |[divM| n = 0 on ri(T,t), 

so that 

Vm ■ Vw dy = [curl u ■ curl w + div u div w] dy . 


(75) 


Comparing dUD and dza, we have that for all test function w e H^{ri{Q,t))j 

[curlM • curl w-f divu divw] dy = • curlrc dy . 


(76) 


We now chose the test function w to have the potential form 

w = Vip , 

for some function periodic function tp e H^(rj(ri,t))- Then, 

curl w = 0 , 


^Note that t) is the image of the 2£-periodic box, and hence functions defined on rjiVl, t) are periodic. 
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and (1761) reduces to 


J div u Atp dx = 0 . 


(77) 


Since u e H^{T}{n,t)) and is periodic, there exists a periodic function tpo e such that 


div u = AijjQ in t). 
Letting ip = ipQ in GZl) then shows that 


(78) 


0 = 


f (divu)2 


dx , 


and thus 


div u = 0 . 

This being true for all time t e [0,T], since t}(x, 0) = x, we then infer that 

detVry = 1. 

Using ((80)) and (IT^ in (|76l) . we see that for all w e 


0 = (£—curl m) • curl w dy , 


and from (IMl) we see that C(v)(x, t) = 0 for all x e fl~, since uJq = 0. Next, we note that 
d,C(v) = = p^A^dvy = Vud) ■ C(v ), 

aXk OXr OXk OXr 

where Vu{r]) denotes Vu o r]. Hence, since C{v) = £ o 77 , it follows that € satisfies 

— Vu • £ = 0 in ry(f2''",t), 


(79) 

(80) 

(81) 


(82) 


and €{y,t) = 0 for all y e ,t)- Since £ G iJ*^( 77 (H+,t)), fc > 2, we take the divergence of 
equation (15^ and find that 


div €t + Vu div £ — Vg; div u + (u^,j d,i —u^i £*,j ) = 0 
From (17^ and the symmetry of the last two terms, we conclude that 

div £( + Vu div £ = 0 , 

and thus 

div £( 77 ( 0 ;, t),t) = div £(a;, 0). 

Since £(0) = curluo we then have from (IMl) that 

div £( 77 ( 0 ;, t), t) = 0 . 

From (l66l) and (l80l) 

f dy,t)dy = 0 

,t) 


(83) 


(84) 

(85) 
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We note that £(-,t) e t)). Next, we define the periodic vector-field V' 6 as the 

solution, modulo constants, of 


—=£ in ? 7 (f 2 +,t), 

—Aip~=0 in ,t), 

with the continuity conditions, which follow from the fact that V'0 e H^{r]{n,t)), 

[V'l = 0 and = 0 on 77(r,<). (86) 

Theorem |4] shows that ip e ,t)), fc > 2. Moreover, from (I84L div-^ is harmonic in both 

r](ri^,t) and r]{n~,t) and is a periodic function; furthermore, |[Vdiv)/; • n| = 0 on r(t), for 

V div Ip- ■ n = curl(curl ip-) ■ n + Aip- ■ n 
= curl (curl ■)/)-) ■ n + ■ n 

= curl(curl ■)/;-) • n 

=Vri(curl■)/’“) • T 2 - Vr 2 (curl^-) • n , 
where we have used £ • n = 0 on r(t) in the third equality, so that 

|[V div^ • n] = [Vri curl'i/j • T 2 — Vra curl'0 • ri| . (87) 


Using (1551) . 
so that 


[curl^] = 0 on? 7 (r,t), 

|[Vi-£, curl^] = 0 on ? 7 (r, t ), 


and from (1571) . 


[V divi/' u] = 0 on7y(r,t). 


( 88 ) 


We now set U-(t) = Using (1551) and the fact that div^ e n iJ*+^(r2-(t)), fc > 2, 

is harmonic in and is a periodic function, we find that 


0 = A div Ip div ip dy + 
Jn+{t) 



A div Ip div ip dy 


r \V div ip\'^dy — f \V div ip\'^dy + 

Jn+{t) Jn-(t) 

r \V div ip\'^dy — f \V div ip\'^dy 

jn+(t) Jn~(t) 


r |[V div div-0 (i5'(t) 
Jr(t) 


which shows that div ip(-,t) is a constant. 

Therefore, 

Alp = — curl(curl 7 />), (89) 

so that — curl(curl■)/)) = £. Substituting this into (I5T]) . we see that for all test functions w e 

0= (—curl(curl'0) — curl m) • curl w dj/. (90) 

Next, we set w = —cuilip -f u in (I90L which satisfies the condition of being a test function, and 
obtain that 

0 = I curl(curl'!/'+ m)P da;, 

J ri{Q,t) 
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and thus 


curlw^ = — curl(curl'!/'^) = £ in , 


and 


curl u = — curl(curl ip ) = 0 in (t) . 
Thanks to (15^ . we have that in 

curl ut + V„ curl u — Vu ■ curl it = 0 , 


which is the same as 


curl {ut + Vuu) = 0 , 


from which we infer the existence of a pressure function p such that 


Ut + Vuit + Vp = 0 . 

Therefore, u is solution of the incompressible Euler equations ©, as we have already proven that 
div It = 0. 

It remains only to show that u{x,0) = uq(x). To this end, we notice that from (I64[) . 

C{v){x, 0) = curlito(a;), 


and thus 


curlit(-,0) = curlitQ, 


which coupled with the fact that divit(-,0) = 0 = divito and the periodicity of it, provides us with 


it(-,0) = ito + c. 


where c is a constant vector. From (l68l) . 


J u{x,0) dx = 0, 

n 

which coupled with 

Uo{x) dx = 0 , 

Jn 

then shows that c = 0 , so that it(-, 0 ) = uo, which completes our proof that it is solution of the vortex 
patch problem on [ 0 ,r], with the desired regularity properties. In particular, by (1541) and (l7T|l . we 
see that ly e ‘^°([0, T]; and hence by the trace theorem, 77 e ‘^°([0, T]; Since 

the vortex patch boundary r{t) = r]{T,t) for each t e [0,T], we see that r(t) is of Sobolev class 
j^k+ 1 / 2 ^ To explain why T(t) is indeed r](T,t), we use the identity curl it 07 ] = Viy • wq, where we 
recall that ojq = curlitg and satisfies (USD. Next, we choose a local coordinate system at a point 
a; e r, such that n{x,0) = 63 and the two tangent vectors are ti = ei and T 2 = 62 - By conditions 
(fT^ .hb we can write Wjj' = X;c(=i ^0 ' SaGa- This means that curlii+(ry(a;, t),t) = rj,a ujq ■ Ca, and as 
we have shown already, curlit+(ry(a;, t), t) • n{ri(x,t),t) = uJq ■ ear],a • ^ Since for a = 1, 2, 

7 y,Q, is a tangent vector to 7 y(r, t) at the point r](x, t) and hence continuous, then 

[curlit] or] = T],a [wcJ" • Ca] . 

This shows that the set r(t), on which curlit(-,t) has a jump discontinuity, is propagated by the 
Lagrangian flow map r](-,t). 

Uniqueness of solutions has been shown by Gamblin & Saint Raymond [8]. 
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6 Elliptic Regularity 


6.1 A two-phase elliptic problem 

For k ^ 2, let c K" denote an open, bounded i7^'''^-doniain which is diffeomorphic to the unit 
ball B = {x : |a;| < 1}. We set T := which is then an it^+^/^-class closed surface. Let 

n denote a periodic box [—£,£]" in with opposite sides identified, and with C sufficiently large 
so that is properly contained in fl. Functions defined on 12 are 2£-periodic in each of the n 
coordinate directions, i.e., 


u{x + 2£ei) = u{x) 'ix e R", f = 1,..., n, 


were Ci denotes the usual Cartesian basis. We set 12 = 12/12+ . 

We establish elliptic regularity for the following two-phase vector-valued elliptic problem: 




du+ \ 


dxj V - dxk ! 


FmI =0 




= 9 
u_ is periodic 


where u+ = ,«!/) and /+ = (/\• • • g = {g^ 

and are two-tensors which satisfy the positivity condition 


in 12-, (91a) 

on r, (91b) 

on r (91c) 

on 512 (91d) 

, g") are vector-valued functions, 


VCeR" (92) 

for some A > 0. We use the notation |[io] = for vector fields w on F, and we let N 

denote the outward unit normal to 512+. The system (1911) has a unique solution in iL^(12) when we 
additionally assume that u(x)dx = 0. 

Let V = iL^(12), the space of functions on [—£, £]” which are 2/1-periodic. Let u = 

Z + /_ln-. The variational (or weak) form of (|M1) is given by 


J 

Jn 




Jni 5xfc dxj 
where we use the following integral notation: 


f fipdx+{ 

Jn± Jr 


dx = ipdx + \ g (pdS y If e V, 


(93) 


and 


f f , f ,fe5wh5¥>*^ 

a-^ -T — — dx = I ai —— -— dx + \ a_ — --— dx 

Jn± Sxk cxj Jq+ dxk cxj Jq- dxk oxj 


f(pdx= \ f+ipdx+\ f^ipdx. 
Jq± Jq+ Jn- 


The regularity theory for solutions u of (IM)) is classical when the coefficient matrix is in 
and can be summarized by the following 

Theorem 3. Suppose that for someke N, e '^^(12+) satisfies (I92p . Then for all e 
and g e the solution u to (1911) is in JL*'+^(12-), and satisfies 


ll^+llir'‘+Hn±) < C'|^||/+||/i-k-i(Q±) -I- ||g||/^k-o.5(r)j 
for some constant C depending on ||a+||‘;j'k(n±). 


(94) 
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We use the following notation for norms: 

ll(■) + llH'‘+l(f^±) = ll(•)+ll^^‘'+l(f2+) + ll(•)-llH'‘+Hf^~) • 

We shall need the corresponding result for the case that the coefficient matrix 0 +' has only 
Sobolev-class regularity: 

Theorem 4. Suppose that for some integer k > ^ and 1 ^ ^ k, e satisfies (1921) . 

Then if f e and g e the weak solution u+ to (1911) is in and satisfies 

\\'>J'±\\H‘+^n±) 

H‘-^{n±) + \\g\\H‘~o-^{r) + ’P(ll«+llff>'(n±)) (ll/+llL2(n±) + ll9llff-o-5(r)^ j : (95) 

where V is a polynomial function and the constant C depends on fl-. 

We are using the notation 

’^(ll®+llH'"(a±)) (ll/+llL 2 (n±) + ||5||//-o.5(r)^ 

= (ll/-i-llL2(n+) + Il9llff-0'=(r)) + 'P{\\o^-\\h'^{q-)) (jif -\\L^(n-) + Il9ll_f/-0'5(r)^ ■ 

Proof. Let E- : iL*'+^(r2-) —> 7L^+^(]R") denote a Sobolev extension operator, and let 0 +^ = 
77 ^ * (E-a+) and = ? 7 e * (E-f). Let {Um}m=i denote an open cover of which intersects the 
interface E, and let {9m}m=i denote a collection of charts such that 

1. 6 m ■ B{0,rm) him IS an iL^+^-diffeomorphism, 

2. det(V0m) > 0, 

3. 6m- B^ = R(0, rm) n {xn = 0} n L, 

4. 6m- B+= B{0,rm) n {j/„ > 0} ^ n , 

5. 6m-B-= B{0,rm) n {j/„ < 0} ^ n , 

6. ||V6» 

771 Id|| 

L’^{B{0,rm)) 1- 

Let 0 ^ Cm < 1 in ^fP{Um) denote a partition of unity subordinate to the open covering Um\ that 
is, 

K 

^ ' Cm ~ 1 and spt(Cm) — him nr. 

771 = 0 

Einally, let denote a smooth regularization of g defined by 

K 

^ g) O 6my\o 6m - 

771=1 

It follows that for e « 1 sufficiently small, 

^ VCeR'.xGff. (96) 

Hence, by Theorem |3l the solution to the variational problem 

f al'"^^—^^dx= r f^ipdxP f g^ipdS VvJeV, 

Jn± oxk dxj Jq+ Jr 
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satisfies u% e for all /c > 1; in particular, the vector fields u% are smooth. We next establish 

an e-independent upper bound for ||||^«+i(f 2 ±p 

Step 1: Regularity in horizontal directions near E. We fix m G {1, and set 

U+ = U%oem, F = f^O dm, G = g^odm, ^ = Cm O 9m , and ^ = (fodm- 

With A = we define 6 ’’® = o9m)AlAj. Then, since || — Id 1 ^ 00 ( 5 +^ « 1, the matrix 

b is positive-definite: 


b^^CrCs = (a^'^ o 9m)AlA^^CrCs > ve e R" . 

By the change-of-variables formula, the variational formulation is written as 

.dx= ( F<^dx+ ( G<^>dS y^eH^{Bm) 
[ Ijrs^ELiAi dx 

JBm dxr dxs 


(97) 


r 

dW d$* 


-7 - 7 — dx = 

Jb± 

dXr dXs 

* j 

— r ,5’"®^^+ 

^ tT<-T 

^ CXs 

JBi; dXr 

^n -1 52 

-ia = l dxf 

r denoting the 


<i>=(-l)^[eA^(ei7)], 


so that 


L 


— dx ^ Cn|/_,_||//f-i(o+) -I- ||/_||ij«-i(n-) + ||5||_f/«-o.5(r) 

i (J Xr (J X a L 


]\\d\CU±) 




(98) 


We focus now on the left-hand side of (1551) . We let d = (di, • • •, dn_i) denote the horizontal 
gradient, and write 

d Ed W ^ ^ dx ■■■dx dx ■■■dx ’ 


d^-Vd^+V =!]••• 2] 


n-l n-1 


dV 


d^AoIE 


ai=l ae-i = l dXa^ ■ ■ ■ dXa^^^ 


and so forth. Then, 

..dC/*d$* 


f d'[b-^(CU),r]d'(CG),sdx- [ d'[b^^UC,r]d'(CU),, dx 

JB^ OXs -ISm -ISm 

+ r ~d'~\b-^U,rC.s]'d'^\mdx. (99) 

■JB^ 

For the first term on the right-hand side of (I99L we make use of (IWl) and Young’s inequality to 
conclude that 

r ~d\¥^{CU),rYd\CU),sdx= \ ¥^t{CU),rt{CU),sdx+ \ [{d\bmcu),r]d\cu),sdx 
JBi, Jb± 
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Then, Corollary |6] with e = 1/8 shows that 

By Lemmaini for 0^£^k+l, /+e and g+ e H^{Q+), and for a generic C, 

II/+ 9±\\H‘(n) < C'l|/+l|ffn'“{k,i>(n±Jg+||//^(n±) 

( 101 ) 

V/+ e e H^n^). 


For the second and third terms on the right-hand side of (1^ , we use the inequality (IIOII) , and find 
that 


f d\¥^u^,r]d\cu),s dx + I d' \b^^U,ri,s]d'^"mda 

^ Cs ||a+ ||l,k(o) II u %-f ^||5 V(Ct/±)). 

Choosing 5 > 0 sufficiently small in fllOOp and f|102p . we conclude that 

II r 1 

ll^d V[/+||^2(B±) ^ ^ ll/+llff«-i(n±) + Il5llff«-0'5(r) + l|a||ffi'(o)||tty|^f+7^^^J . 

Step 2: Regularity in the vertical direction near T. We write (I^Tk l as 

-^{¥^U+,s),r=^F+ in R±. 


(102) 


(103) 


(104) 


We analyze (11041) in the -f-phase and drop the -I—subscript for notational clarity. With !/,„ denoting 
dUIdxn, we have that 

=^[F- b--,n C/,„ - C. - sr ] in B+ . (105) 

We analyze the terms on the right-hand side of (jl05l) . For any integer j such that 0 ^ j ^ — 1, 

< C[||/||„,. 

Moreover, since £ ^ k, by Lemma[5]with e = 1/8, 


||d' ' "v^(e6'^",„£/,„)||^2(B+)+ E 

(r,s)^(n,n) 




i-1 


^ C y l|V^-’-aD’'+iu^|U2(aF) ^ C y ||V^+i-^ai^’-n^| 

r=0 r=l 

Finally, by Corollary [5] with e = 1/8, 


L2(0 + ) ^ Ce||a||/fk(Q+/|tt ■ 






(r,s)75:(n,n) 

^ C'£||a||ffk(o)||u'=||^,_^7^^^^. 


Therefore, for 0 ^ j ^ £ — 1, letting d^ ^ act on ()105F 




(106) 


(r,s)7^:(n,n) 
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for a function G(£j) satisfying 

Now we argue by induction on 0 ^ j ^ ^ — 1. By (EH), 6™ ^ -I so that when j = 0, the inequalities 
(11031) and (I106|) show that 

1^2(5+) ^ II G(^j)||^2(^+) + 2 II& IIl®(B+)||?^ IIl2(b+) 


(r,s)^(n,n) 


which, combined with (I103L provides the estimate 




u 


(n^ 


Repeating this process for j = 1, • • • , £ and including the analysis in the —phase, we conclude that 


^ C \\f ±\\H‘-^{n±) + ll5llff^-o-5(r) + l|a+||_f/k(Q±)||My|^f+7^^^J . 


(107) 


Step 3: Completing the regularity theory. Let y+ > 0 be in so that spt(x+)ccifl-. 

Repeating the computations above, we find that 

||x+V^+^uy|i2(Q±) ^ c|^||/+||ff<!-i(o±) + ||a+||ijk(n±)||u+||^j+7^^^J . (108) 

The inequalities (I107p and (11081) establishes the inequality 

ll^^+llff«+i(n±) < G|^||/+||/i-f-i(n±) + ||ff||/i-^-o.5(r) + l|a+llHk(Q±)||u+l|^j+7^^^J . (109) 

Since 

II'“+II^<!+|(q+) ^ G||M+||^4)(n±)ll'“+lllfi(n+) ’ 

Young’s inequality shows that 
ll^±llff'!+i(a±) 

< G5|^[|/+||^f-i(Q±) + ||5||/i-^-o.5(r) + 'P(||a+||^k(Q±))||it+||//i(Q±)j + (5 ||m+||^^+i(o±) 

for some polynomial function V. Finally, the inequality (IMl) is established by choosing i5 > 0 
sufficiently small, letting e —> 0, and using the a priori iJ^-estimate. □ 


A Some basic inequalities 

Lemma 5. For k > ^ and 0 ^ ^ k, let O R" be a bounded smooth domain. Then for all 

e G (O, i), there exists a constant depending on e such that for all f e H^(0) and g e H^~^{0), 

i 

X ||VVV'-^5||i2(o) < G,||/||^k(o)||5ll//-qo) • (110) 

1=1 
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Proof. We estimate ^<7 for j = 1, • • • , £ as follows: 

Step 1. If 1 ^ j ^ by the Sobolev inequalities 


we find that 


(if 0 < e< 1), 


L„- 2 U -0 (O) 


l|VVV^-5lk2(o) < 

^ C'e||/||^§+ej.Q^||5||//i-qO) • 


Step 2. If j = £, by the Sobolev inequality 


we find that 


I|w^I|l<=°(0)< C'<;||u;||^n+e|-Q^ , 


Step 3. If ^ < j < £ (this happens only when ^ ^ k), we consider the following two sub-cases: 

Case ^ n: Similar to the previous case, by the Sobolev inequalities 

C'||w||^§-f+j^Q^, 

and hence, we obtain that 


||w|| _ 2n ^ C\\M\H^-i(0) and ||w| 


(O) 


VVV^-^5lU2(o) < < C||/||^.(0)ll5ll^S 


L^(0) 


"H2(0) ' 


Case B: n < £ ^k: If j > k — - , by the Sobolev inequalities 


||w|| _ 2 n_^ ^ C'l|w||_H-k-3(o) and ll'icl 


-L’^(O) 




H7-'‘+2(0) 


we obtain that 


< C'll/llff>2(o)ll5ll ui -k-b^ 




l|VVV^-5lU2(0) ^ 

Now suppose that ^ ^ ^ ^ ^ ^ ■ Note that if 0 < e < i, 

+ ^ C'e||w^llwJ.«>(0) ^ C'e||^llff>2(0) > 

Therefore, by the Gagliardo-Nirenberg-Sobolev interpolation inequality, we obtain that 
l|V^fV^-^g||i2(o) < ll/ll wj'”(0)I|5||h'^-j(0) 
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for some aj G (0,1); hence, by Young’s inequality, 


llWV'-^gll L2(o) < <^£[11/11 jjt+=(o)ll5llH^-dO) + ll/llffi‘( 0 )ll 5 ll^t-k+^( 0 )] ■ 


Summing over £, we conclude that for 0 < e < -, 







otherwise 


Estimate (IllOl) is then obtained from the fact that for all e e (O, i), 


^ + e ^ k and ^ + e ^ ^ — e if (in addition) ^ ^ • 


□ 


Corollary 6. For any m G {1,..., K}, and for F G and G = with 0 < e < 1/4 

and 1 ^ ^ k, 

C4F\\ \\G\\(Bf,) ’ (111) 

where {d^F}G = d^FG) - Fd^G. 
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